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ABSTRACT 

By means of cosmological N— body + Hydrodynamics simulations of galaxies in the context of the A Cold 
Dark Matter (ACDM) scenario we explore the specific star formation rates (SSFR=SFR/M S , M s is the stellar 
mass) and stellar mass fractions (F s = M s /Mh, Mh is the halo mass) for sub-M* field galaxies at different 
redshifts (0 ^ z < 1.5). Distinct low-mass halos (2.5 £ M h /10 W M Q < 50 at z = 0) were selected for the 
high-resolution re-simulations. The Hydrodynamics Adaptive Refinement Tree (ART) code was used and some 
variations of the sub-grid parameters were explored. Most simulated galaxies, specially those with the highest 
resolutions, have significant disk components and their structural and dynamical properties are in reasonable 
agreement with observations of sub-M* field galaxies. However, the SSFRs are 5-10 times smaller than the 
averages of several (compiled and homogenized here) observational determinations for field blue/star-forming 
galaxies at z < 0.3 (at low masses, most observed field galaxies are actually blue/star-forming). This inconsis- 
tency seems to remain even at z ~ 1 — 1.5, although it is less drastic. The F s of simulated galaxies increases 
with Mh as semi-empirical inferences show. However, the values of F s at z « are ~ 5 — 10 times larger 
in the simulations than in the inferences; these differences increases probably to larger factors at z ~ 1 — 1.5. 
The inconsistencies reported here imply that simulated low-mass galaxies (0.2 ;$ M S /10 9 M Q ;$ 30 at z = 0) 
assembled their stellar masses much earlier than observations suggest. Our results confirm the predictions 
fou nd by means of ACDM-based models of disk galaxy formation and evolution for isolated low-mass galax- 
ies dFirmani & Avila-Reesell20Tol) . and highlight that our understanding and implementation of astrophysics 
into simulations and models are still lacking vital ingredients. 
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1. INTRODUCTION 

Based on the recent observational achievements, an empir- 
ical picture of stellar mass (M s ) assembly history of galaxies 
as a function of mass is emerging. Probably, the most remark- 
able aspect in this picture is the so-called 'dow nsizing' phe- 
nomen on. The term 'downsizing' was coined bv lCowie et al.l 
d 19961) to describe the rapid decline with cosmic time of the 
maximum rest-frame A'— band luminosity of galaxies under- 
going active star formation (SF). Since then, as new observa- 
tional inferences appeared, this term has been used to describe 
a number of trends of the galaxy population as a function of 
mass, most of them related actually to different astrophysi- 
ca l phenomena and gal axy evolutionary stages as discussed 
in iFontanot etall(l2009l see also Neistein et al. 2006; Santini 
et al. 2009) . 

From the most general point of view, the many downsiz- 
ing manifestations can be separated into those that refer to 
the evolution (i) of massive galaxies, which today are on av- 
erage red and passive (quenched SF), and (ii) of less mas- 
sive galaxies, which today are on average blue and star form- 
ing. This separation seems to have an astrophysical root; 
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that is, it is not based on a merely methodo logical division. 
For example, iFirmani & Avila-Reesd (1201 Ol hereafter FA 10) 
have inferred the average M s growth tracks of galaxies as a 
function of mass (the 'galaxian hybrid evolutionary tracks', 
GHETs) by means of a semi-empirical approach, and showed 
that at each epoch there is a characteristic mass that separates 
galaxies into two populations. Galaxies more massive than 
M s (z = 0)w 3 x 10 10 M Q are on average passive (their SS- 
FRs have dramatically diminished, where SSFR=SFR/M S is 
the specific SF rate), besides the more massive is the galaxy, 
the earlier it transited from the active (blue, star-forming) to 
the passive (red, quenched) population ('population downsiz- 
ing'). Galaxies less massive than M s (z = 0)~ 3 x 10 10 
Mq are on average yet blue and active, and the less massive 
the galaxy, the higher its SSFR and/or the faster its late M s 
growth ('downsizing in SSFR'). We recall that these are just 
average trends. In fact, there are other factors besides mass 
that intervene in the stellar mass build-up of galaxies, e.g., 
the large-scale e nvironment in which galaxies evolve (e.g., 
lPenget alJl20Tol and more references therein) and whether 
they are central or sat ellites objects (for a recent review see 
IWeinmann et al.ll2011l) . 

A large amount of direct look-back time observations sup- 
port the mentioned two downsizing phenomena: 

• For high-mass galaxies, those with M s ^ M*, where 
M* is the characteristic mass of the stellar mass func- 
tion (w lO 10 8 M Q at z ~ 0), observations show the 
existence of a decreasing with cosmic time characteris- 
tic mass at which the SFR is dramatically quenched or 
at which the stellar mass functions of early- and late- 
type galaxies cross, evidencing this a mass-dependent 
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migration from blue to red population (e.g., Bell et al. 
2003,2007; Bundy et al. 2006; Hopkins et al. 2007; 
Drory & Alvarez 2008; Vergani et al. 2008; Pozzetti et 
al. 2010). 

• For low-mass galaxies (M s << Al*, they are 
mostly late-type, star-forming systems), at least up 
to z ~ 1 — 2, and by using different SFR tracers 
and methods to estimate M s , it was found that their 
SSFRs are surprisingly high even at z ~ and, on 
average, the lower the mass, the h igher th e SSFR (e.g.. 



Baldrv et al.l l2004t lBauer etal.1 120051: iFeulner et al l 
20051; IZheng et al. I 120071; iNoeske et al. I l2007bT 



Bell et al. I l20Q7t lElbaz etafl 2007; Sa lim et alj 120071; 



Chen et al.| 120091; iDamen et al.1 l2009allbl; ISantini et all 



20091: lOliveretalJ 12010b iRodighiero et al.l l2010t 



Karim et al.ll201 lh 



1.1. Confronting the empirical picture to theoretical 
predictions 

The current paradigm of galaxy formation a nd evolution is 
based on the hierarchical c lustering scenario (IWhite & Reesl 
ll978HWhite & FrenklfT99Tl) within the context of the A Cold 
Dark Matter (ACDM) cosmology. According to this scenario, 
galaxies form and evolve in the centres of hierarchically grow- 
ing ACDM halos. Though at first glance contradictory, it 
was shown that the 'population (or archaelogical) downsiz- 
ing' related to massive galaxies has, at least partially, its nat- 
ural roots in the hierarchical clustering proces s of the dark- 
matte r halos and their progenitors distribution dNeistein et alJ 
120061 see also Guo & White 2008; Li, Mo & Gao 2008; Keres 
et al. 2009). Besides, red massive galaxies are expected to 
have been formed in early collapsed massive (associated to 
high peak and clustered) halos that afterwards become part 
of groups and clusters of galaxies, leaving truncated therefore 
the mass growth of the galaxies associated to these halos. On 
the other hand, massive galaxies typically hosted in the past 
active galactic nuclei (AGN). The strong feedback of the AGN 
may help to stop gas accretion, truncating further the galaxy 
stellar growth and giving ris e to shorter formation time-scales 
for more massive galaxies dBower et al.l 120061: ICroton et al.l 
120061 iDe Lucia et al.ll2006h . 

All the factors mentioned above work in the direction of re- 
producing the downsizing manifestation s of massive galaxies 
within the hierarchical ACDM scenario dFontanot et al.l2009h 
-though several questions remain yet open. 

What about t he downsizing related to low-mass galaxies? 
As discussed in iFirmani et al.l (1201 Ol) and FA 10, disk galaxy 
evolutionary models in the context of the hierarchical ACDM 
scenario seem to face difficulties in reproducing both the high 
empirically determined values of SSFR and the SSFR down- 
sizing trend of low-mass galaxies (M s £3 x 1O 1O M << 
M*), which are mainly blue/star-forming systems of disk-like 
morphology. 

The gas infall rate onto the (model) disks is primar- 
ily driv en by the cosmological halo mass aggregation 



history (|Fir mani & Avila-Reesel l2000t: Ivan den Boschl 1 2000 ; 
Stringer & Benson! 120071: IDutton & van den Boschl 12009 



Firmani et al.ll20l6T: IDutton et al.ll2010ah . which is hierarchi- 
cal, i.e. the less massive the halo, the earlier its fast mass 
aggregation rate phase. The gas infall rate is the main driver 
of disk SFR (other factors related to the process of gas conver- 
sion into stars, local feedback, interstellar medium turbulence, 
etc. introduce minor systematic deviations). The large-scale 



feedback over the gas exerted mainly by SNe explosions, may 
strongly affect the SFR history of low-mass disks mainly be- 
cause it alters the primary gas infall rate. The SN feedback has 
been commonly invoked in semi-analytical models (SAM) to 
reproduce the faint-end flattening of the galaxy luminosity (or 
M s ) function, assuming that this feedba ck produces gas re- 
heating and galactic outflows (see e.g., iBenson et al.1 120031 
and more references therein ). IFirmani et al. I (12010b experi- 
mented with different models of SN-driven galactic outflows 
(e.g., based on energy and momentum conservation) besides 
of the disk gas turbulence input due to SN feedback. They 
have found that the present-day stellar and baryonic mass 
fractions (F s = M s /M h and F b = M b /M h , where M h is the 
total -virial- halo mass and M b = M s + M g , M g is the galaxy 
cold gas mass) of low-mass galaxies can be roughly repro- 
duced when assuming extrem e (probably unrealistic) galactic 
outflow efficiencies (see also IDutton & van den Boschll2009l; 
IDutton et al.l 1201 0"ah . However, in any model with galactic 
outflows it was possible to alter the SFR histories of their 
disk galaxies in order to reproduce the observed (at least up to 
z ~ 1) downsizing in SSFR and the too high SSFRs measured 
f or small central galaxi es at late epochs. 

IFirmani et al.l (1201 Ol) experimented also with the possi- 
bility of late re-accretion of the ejected gas (in previ- 
ous works it was assumed that the feedback-driven out- 
flow s eject gas from t he small halos forever, but see 
e.g..lSpringel et al.ll2001t IDe Lucia et al.ll20"04t iBertone et al] 
120071: lOppenheimer & Davell2008l:IOppenheimeret alj|2010l) . 
For reaso nable schemes of gas re-accretion as a function of 
halo mass. IFirmani et aD(l2010l) found that the SSFR of galax- 
ies increases but it does it in the opposite direction of the 
downsizing trend: the increase is larger for the more massive 
galaxies. 

A related issue of ACDM-based low-mass (disk) galaxy 
models appears in the evolution of the M s -Mh (or F s - 
Mh) relation. The semi-empirical inferences of this rela- 
tion (see e.g IConrov & Wechsler l 120091: iMoster et al.ll2010l 
IWang & Jind2010l;lBehroozi et al.l2010l) show that for a fixed 
value of Mh, F s is not only very small at z = 0, but it be- 
comes smaller at higher z. The ACDM-based disk galaxy 
evolutionary models predict the opposite, i.e. higher values of 
F s at higher z for a fixed value of Alh (FA10). 

The inconsistence between models and observations related 
to the stellar mass bui ld-up of sub-M* disk galaxies found 
in IFirmani et al.l (1201 Ol) and FA 10 is connected with some is- 
sues reported in recent SAMs. In these works, based also 
on the hierarchical ACDM scenario, it was showed that the 
stellar population of small (both central and satellite) galaxies 
(M s ss 10 9 — 10 10 5 M ) is assembled too early, becoming 
these galaxies older, redder, and with lower SSFRs at later 
epochs than the observ e d galaxies in the same mass range 
dSomerville et all 120081: iFontanot et all 120091: ISantini et al.1 
120091; iPasquali et al.ll2010t iLiu et al.l |20 1 Ol) . B y means of a 
disk ga laxy evolutionary model similar to FAlO. lDutton et all 
(1201 Oal) found also that the SSFR of low-mass galaxies is be- 
low the average of that given by observations, specially at 
z ~ 1 — 2, though these authors conclude that their mod- 
els roughly reproduce the main features of the observed SFR 
sequence. 

1 .2. Cosmological simulations of galaxy evolution 

Due to the high non-linearity implied in the problem 
of halo/galaxy formation and evolution, cosmological N- 
body/hydrodynamical simulations offer the fairest way to at- 
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tain a realistic modeling of galaxies. However, this method 
is hampered by the large -currently unreachable- dynamic 
range required to model galaxy formation and evolution in the 
cosmological context, as well as by the complexity of the pro- 
cesses involved, mainly gas thermo-hydrodynamics and SF 
and its feedback on the surrounding medium. The latter pro- 
cesses occur at scales ( ;$ 1 — 50 pc) commonly below the 
accessible resolution in simulations. Therefore, "sub-grid" 
schemes based on physical models for these processes should 
be introduced in the codes. 

In the last years, cosmological simulations have ma- 
tured enough as to produce indiv idual galaxies that at 

z ~ look qu i te rea li stic (e.g.. iGovernato et all 120071 

2010t lNaabetal.1 120071: iMaver et all 120081; IZavalaetal l 



20To| B 

TOOa ICeverino&Klypinl 120091: bibson et all 1 20091: 



Scannapieco et al.l 2008: |Piontek & Steinmetzl 12011 



Colin et all I2oToF ISawala et al.l 1201 It lAgertzetall I20TI 



though this success in most cases is conditioned by the 
assumed sub-grid schemes and parameters, and the resolution 
that limits the simulations to small boxes, where only one or 
a few galaxies are followed. 

What do current cosmological N-body/hydrodynamical 
simulations of low-mass galaxies predit regarding their SSFR 
and F s evol ution? By using the Adaptive Refinement Tree 
code ( ART IKravtsov et al] 1997[) w i th hyd rodynamics in- 
cluded dKravtsovl 120031) . IColm et alJ (1201 Ol hereafter C10) 
have explored different SF schemes and sub-grid parameters 
for one simulated low-mass galaxy (Mh ~ 7 x 10 10 M© at 
z = 0) that develops a significant disk. For some cases, the 
obtained galaxy did not look like realistic. For other cases, 
several global dynamical, structural, and ISM properties of 
observed small galaxies were reproduced. However, even in 
these cases the SSFR and F s at different z are respectively 
much lower and higher than the observational inferences, an 
issue that apparently would also share other simulations pre- 
sented in the literature. 

In this paper, our aim is to measure at different epochs 
the SSFRs and stellar mass fractions (F s ) of sub- M* central 
galaxies covering a large range of masses obtained in state- 
of-the-art N-body/hydrodynamical cosmological simulations, 
and explore whether the mentioned-above issues of ACDM- 
based galaxy evolutionary models (as well as SAMs) are also 
present or not in the simulations. Here we will use the Hydro- 
dynamics ART code with some of the most reliable sub-grid 
schemes/parameters explored in C10 and covering almost two 
orders of magnitude in M s , The effects of resolution over the 
SSFR and F s at different epochs will be also explored. 

The plan of the paper is as follows. The code and the 
SF/feedback prescriptions used for the simulations are de- 
scribed in §2. In this section, the cosmological simulation and 
the different runs varying the SF/feedback parameters are also 
presented. In §3.1 some generalities of the simulated galaxies 
are discussed, while in §3.2 the SSFRs and F s of all the simu- 
lations at z — 0.00, 0.33, 1.00, and 1.50 vs mass are presented 
and compared with several observational inferences. In §§4.1 
we discuss the numerical results from previous works, and 
in §§4.2 the caveats of the observational determinations used 
here are discussed. Our conclusions are presented in §5. 

2. THE SIMULATIONS 
2.1. The Code 

The numerical simulations used in this work were per- 
formed using the Hydrodynamics + N-body ART code 



( IKravtsov et al.lll997t lKravtsovll2003h . Among the physical 
processes included in ART are the cooling of the gas and 
its subsequent conversion into stars, thermal stellar feedback, 
self-consistent advection of metals, a UV heating background 
source. 

The cooling and heating rates incorporate Compton heat- 
ing/cooling, atomic, and molecular Hydrogen and metal-lines 
cool ing, UV heating from a cosmological background radia- 
tion (lHaardt&Madaull99 61 and are tabulated for a tempera- 
ture range of 10 2 < T < 10 9 K and a grid of d ensities, metal- 
licitie s, and redshifts using the CLOUDY code dFerland et al. I 
119981 version 96b4). 

Star formation is modeled as taking place in the coldest 
and densest collapsed regions, defined by T < Tgp and 
n g > rigp, where T and n g are the temperature and num- 
ber density of the gas, respectively, and Tgp and ngp are a 
temperature and density SF threshold. A stellar particle of 
mass m* is placed in a grid cell where these conditions are 
simultaneously satisfied, and this mass is removed from the 
gas mass in the cell. The particle subsequently follows N- 
body dynamics. No other criteria are imposed. In most of 
the simulations presented here, the stellar particle mass, m*, 
is calculated by assuming that a given fraction (SF local effi- 
ciency factor csf) of the cell gas mass, m g , is converted into 
stars; that is, m* = egpm s , where £sf is treated as a free pa- 
rameter. Based on the simulations performed in C10, where 
actually egp was not a fixed parameter but a quantity calcu- 
lated in each cell and timestep by an algorithm dependent on 
other parameters, we have found that when egF acquires val- 
ues around 0.5, the simulation gives reasonable results. This 
value is high enough for the thermal feedback to be efficient 
in regions of dense, cold, star-forming gas and not as large as 
to imply that most of the cell gas is exhausted. 

In C10, besides the "deterministic" SF prescription de- 
scribed above, the authors experimented also with a "stochas- 
tic" or "random" SF prescription in which stellar particles are 
created in a cell with a probability function proportional to the 
gas density. The stochastic prescription allows for the possi- 
bility (with low probability) of forming stars in regions of low 
average density, in which isolated dense clouds are not re- 
solved. Here, all of our simulations, except one, use the "de- 
terministic" SF prescription assuming a density threshold of 
ngp = 6 cm" 3 , which corresponds to a gas column density 
N ~ 10 21 cm -2 (roughly the lower limit of observed giant 
molecular clouds) in a cell of 100 — 150 pc (see C10). 

Since stellar particle masses are much more massive than 
the mass of a star, typically 10 4 - 10 5 M©, once formed, 
each stellar particle is considered as a single stellar popula- 
tion, within which the individual stellar masses are distributed 
according to the Miller & Scalo IMF. Stellar particles eject 
metals and thermal energy through stellar winds and type II 
and la supernovae (SNe) explosions. Each star more massive 
than 8 M is assumed to dump into the ISM, instantaneously, 
2 x 10 51 ergs in the form of thermal energy; 10 51 ergs comes 
from the stellar wind, and the other 10 51 ergs from the SN 
explosion. Moreover, th e star is assumed to eje ct 1.3M© of 
metals. For the assumed [Miller & Scalo I (1 19791) initial mass 
function, IMF, a stellar particle of 10 5 M© produces 749 type 
II SNe. For a more detail ed discuss i on of the p rocesses imple- 
mente d in the code, see iRravtsovl d2003l) and IKravtsov et alJ 
d2005h . 

Stellar particles dump energy in the form of heat to the cells 
in which they are born. Most of this thermal energy, inside 
the cell, is radiated away unless the cooling is turned off tem- 
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porally. This mechanism along with a relatively high value 
of esF allow the gas to expand and move away from the star 
forming region. Thus, to allow for outflows, in the present 
paper we adopt the strategy of turning off the cooling during 
a time t a s in the cells where stellar particles form (see CIO). 
As t g can be linked to the crossing time in the cell at the 
finest grid, we could see this parameter as depending on reso- 
lution in the sense that the higher the resolution, the smaller its 
value. In most simulations presented in this work, we turn off 
the cooling for t Q s = 40, 20, and 10 Myr depending on reso- 
lution; however, we have found that varying this time, within 
this range, for a given simulation, does not affect significantly 
the results. 

2.2. Numerical strategy and the runs 

Most simulations presented here were run in a ACDM cos- 
mology with ft = 0.3, n A = 0.7, n b = 0.045, and h = 0.7 
The C DM power spectrum is taken from lKlvpin & Holtzman I 
( 1 1997b and it is normalized to erg = 0.8, where a% is the rms 
amplitude of mass fluctuations in 8 /i _1 Mpc spheres. Few 
simulations were run using a different cosmological set up 
(models C, F, and H of Table 1). For them, fi = 0.27, 
A = 0.73, 6 = 0.047, and h = 0.7, and the power 
spectrum is that on e used to run the "Bolshoi simulation" 
dKlvpin etal. 11201 Oh with er 8 = 0.82. 

To maximize resolution efficiency, we use the "zoom-in" 
technique. First a low-resolution cosmological simulations 
with only dark matter (DM) particles is run, and then regions 
(DM halos) of interest are picked up to be re-simulated with 
high resolution and with the physics of the gas included. The 
low-resolution simulations were run with 128 3 particles in- 
side a box of 10/i _1 Mpc per side, with the box initially cov- 
ered by a mesh of 128 3 cells (zero-th level cells). At z = 0, 
we search for low-mass halos (2 x 10 10 M/ 1 //i _1 M© ;$ 4 x 
10 11 ) that are not contained within larger halos (distinct ha- 
los), but the selection was not based on their MAHs. 

Out of ten halos, seven do not have companion halos at z = 
with a mass greater than 0.2 the mass of the selected halo 
within a distance 0.5/i _1 Mpc. On the other hand, two of 
them are 220 /i _1 kpc away from each other (halos G and J, 
see Table 1), and finally, there is one halo (D) that is separated 
by 330 /i _1 kpc from a halo with a mass of about 0.5 the mass 
of halo D. 

A Lagrangian region of 2 or 3i? v ir is identified at 
z = 50 and re-sam pled with additional small-scale modes 
dKlvpin et al. Il200l"l) . The virial radius, i? v i r , is defined as the 
radius that encloses a mean density equal to A v j r times the 
mean density of the universe, where A v i r is a quantity that 
depends on Q,q, Q\ and z. For example, for our cosmology 
A vil (z = 0) = 338 and A vir (z = 1) = 203. The number 
of DM particles in the high-resolution region depends on the 
number of DM species and the mass of the halo, but for mod- 
els with four species (high-resolution) this vary from ~ half a 
million (model Ah) to about 4.7 million (model Ih), the least 
massive and the second most massive halos, respectively. The 
corresponding dark matter mass per particle of model galaxies 
(wdm) is given in Table 1 . 

In ART, the initially uniform grid is refined recursively as 
the matter distribution evolves. The criterion chosen for re- 
finement is based on gas or DM densities. The cell is refined 
when the mass in DM particles exceeds 1.3(1 — F bt u)m p or 
the mass in gas is higher than l3.0Ft,u m p , where m p = 
m dm + m g = m dm /(l - F bjU ) and F b)U is the universal 
baryon fraction, F b ,u= 0.15 for the cosmology used here; it 
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Parameters of the simulations 
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"Stochastic SF scheme 

^"Run away" model, without turning off cooling 

is assumed that the mean DM and gas densities in the box 
are the corresponding universal averages. For the simulations 
presented in this paper, using multiple DM particle masses, 
the grid is always unconditionally refined to the third level 
(fourth level), corresponding to an effective grid size of 512 3 
(1024 3 ). On the other hand, the maximum allowed refinement 
level was set to 9 and 10, for low and high resolution simula- 
tions, respectively. This implies spatial sizes of the finest grid 
cells of l max = 218 and 109 comoving pc, respectively. 

We have re-simulated, as described, the evolution of ten 
individual galaxy/halo systems covering a total halo mass 
(baryon + dark matter) range from M/, « 2 x 10 10 M Q to 
Alh ~ 4 x 10 11 M Q at z = 0. For most of these systems, we 
have run two or more simulations by varying the resolution, 
the SF local efficiency parameter esF, the time during which 
cooling time (after SF events) is kept off, < D ff, and, in one 
case, the SF scheme. Table [T] summarizes the characteristics 
of all the simulations studied here. 

The first column gives the name of the run, denoted by 
a capital letter that refers to a particular galaxy /halo system 
(from A to J for ten systems), followed by a lowercase let- 
ter that is labeled "I" or "h" meaning low (218 pc) and high 
(109 pc) resolution, respectively, and finally by a number that 
is used to differentiate cases with different combinations of 
the parameters esf and t Q ff (indicated in the next columns). 
There are three simulations marked with the symbol "*". In 
simulations with "h" but without "*", cells were refined if 
they satisfy the above mentioned condition but with m p given 
by its lower resolution (less species) value, which is 8 times 
larger. Thus, the refinement condition for these simulations 
poorly obeys the cell scaling and the "h" simulations end up 
with a number of resolution elements close to that obtained 
actually for "I" simulations. Model galaxies with "h" and 
"*" do satisfy the scaling: cells are refined according to the 
above recipe with m p given by its higher resolution value. In 
other words, simulations with "h" and "*" are of high resolu- 



Cosmological simulations of low-mass central galaxies 



5 



o 

-0.2 

o 

-0.4 

sf 

60 -0.6 
o 

-o.a 
i 



-0.2 

o 

? -0.4 

SO —0.6 
o 

-0.8 

-1 


-0.2 
\ -0.4 
g 1 -0.6 
0.8 

1 



_ 


1 II 1 1 1 

1.9<M h0 <3.0 

\ 

X 


- — 

"a v \ 

- A \ _x 

- — c > V \ - 

EPS(3.2) \ 




3.0<M h ,o<8.0 


7 D 

- E 

7 F 

EPS(6.[ 


\ ^^\~x. 

\ V 

i) \ V 

1 ""' 




8.0<M h0 <40.0 


\ \ ^x, 

- \ \ _ 

\ \ 
x \ x 

\ \ \ x. 
7 G \S s \ - 

I \\ \ \ 

1 — J V \ \ \ - 

EPS(25.1) \ \ 

1 1 \ l\ X l\ 1 1 1 



3 4 
1 + z 



6 7 8 9 10 



FIG . 1 . — Mass aggregation histories of nine of our simulated halos normal- 
ized to the mass attained at z = 0, q. In each panel, the MAHs of three 
halos (identified with the uppercase letter) in the mass range shown in the 
upper-right parts of the panels are plotted (Afh.o is m unities of 10 10 Mq). 
The dashed curves correspond to the average MAHs (20,000 realizations) of 
halos of Mho as indicated in the parenthesis. These MAHs were calculated 
using the EPS formalism by Firmani & Avila-Reese (2000). The horizontal 
dotted lines correspond to half - 



tion in the sense of the number of DM particles and in terms 
of number of grid cells. As can be anticipated, the modeling 
of galaxies with "h" and "*" consumes a lot of CPU and wall 
time and thus, although desirable, not all models could be run 
with this resolution. In columns (2), (3), and (4) the total halo 
virial mass Mh at z = 0, the DM particle mass resolution, 
rridm, and the size of the finest grid cell, £ max , are reported. 
The values used for the parameters esf and t a s are shown in 
columns (5) and (6). 

Galaxies G and J are actually from one simulation, corre- 
sponding to the box used in CIO. This simulation was run to 
get a maximum resolution of £ max = 218 pc (three DM par- 
ticle species). Cases in which the egp column is written " — " 
refer to the simulations where egp was not fixed but was cal- 
culated at each cell and timestep by an algorithm introducing 
other parameters (see §§2.1 and CIO). The system E is the 
one that covers the largest variation of cases and parameters. 

Galaxy H was simulated with the stochastic approach for 
SF, csf not fixed a priori, and without turning off the cool - 
ing. In this simulation, following lCeverino & Klvpiril ( 120091) . 
"runaway stars" are modeled. These are massive stars with 
high peculiar velocities as a result of the SN explosion in a 
close binary or due to two-body encounters in stellar clus- 



ter. As a result, runaway stars can move far away from the 
SF regions, being then able to inject energy into the ISM in 
regions of low-density gas, where the cooling time is large 
so the feedback efficiency is higher. The runaway stars are 
modeled through the addition of a random velocity, drawn 
from an exponential distribution with a characteristic scale of 
17 km s~ , to the inherited gas velocity in 30% of the stellar 
particles. 

2.3. Halo mass aggregation histories 

The halo mass aggregation histories (MAHs) of the sim- 
ulations presented in Table 1 , except for run H , are plotted 
in Fig. Q] The SF histories of galaxies are expected to fol- 
low in a first approximation the halo MAHs. As seen, the 
MAHs of the different simulated distinct halos are quite di- 
verse, and several of them oscillate around the average ACDM 
halo MAHs expected for these masses (dashed curves cal- 
culated with the Extended Press-Schechter -EPS- formalism 
of Firmani & Avila-Reese 2000, which gives results in good 
agreement with the fits to the outcomes of the Millenium and 
Millenium II simulations (Fakhouri, Ma & Boylan-Kolchin 
2010) 6 . Most of our simulated halos show late (since z ~ 0.6) 
MAHs close to the corresponding averages. 

Halos B, F, and J suffer a late major merger at z « 1.0, 
1.1, and 0.7, respectively. Actually, halo F suffers two major 
mergers that occur around z = 1; this explains the big jump 
in mass seen in the MAH of this halo. Halos A and E have 
MAHs that imply very early mass assembly. Halos C, F, I 
and J have MAHs implying mass assemblies later than the 
corresponding averages. In particular, is remarkable the de- 
layed MAH of halo C. The increase in mass seen at z k, 0.33 
in this halo is due to the almost simultaneous accretion of 
three medium-sized halos each of mass of about one tenth of 
the mass of halo C, two of them coming in a pair. 

It is important to remark that all the simulated halos are 
distinct at z = 0, and the galaxies studied here are the central 
ones inside these halos. 

3. RESULTS AND COMPARISON WITH 
OBSERVATIONS 

Most of the simulated galaxies studied here show dynam- 
ical, structural, and ISM properties in reasonable agreement 
with observations of sub-M* central galaxies. The code 
SF/feedback schemes and their parameters (see §2.1) were ex- 
tensively explored in C10. Based on the analysis results from 
that work, here we limit our experiments to only a few varia- 
tions of some of the parameters and to simulations with higher 
resolutions (see Table[T]i. Our goal is to explore the SFR, M s 
and Mh evolution of low-mass isolated galaxies in the numer- 
ical simulations and to compare them to current observational 
inferences of central sub-M* galaxies at different redshifts. 

3.1. General properties 

The main properties at z = of all runs are reported in 
Table |2] The simulation name and halo virial mass Mh are 
given in columns (1) and (2). Columns (3) and (4) present the 
galaxy stellar (disk + spheroid) mass and the (cold) galaxy 
gas mass fraction, M s and f g =M g l{M s + M g ), respectively, 
where M g is the mass contained in gas particles with T < 10 4 

6 While the EPS average MAHs are for pure DM halos, the MAHs from 
the simulations are for the DM+baryons systems. The inclusion of baryonic 
physics actually does not affect significantly the halo MAHs. 
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Physical properties at z = of simulated galaxies 
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K. Both M s and M g are defined as the stellar and gas mass 
contained within 5R e , where R e is in turn defined as the ra- 
dius where half of the galaxy stellar mass is contained (col- 
umn 5). This latter mass is found inside the radius defined 
as the minimum between the tidal and 0.5 the halo virial ra- 
dius. This definition considers the stellar mass associated with 
a possible very extended stellar halo; satellite galaxies would 
be also considered, but this makes no difference as far as R e 
and related properties are concerned because their contribu- 
tion to the total mass is ^ 5%. The maximum circular ve- 
locity, V max , is given in column (6); the circular velocity is 
computed as V c = y / GM(R)/R, where M{R) is the total 
(dark, stellar, and gas) mass. Finally, columns (7) and (8) 
give the galaxy stellar and baryon mass fractions, defined as 
F a =MJM h and F b =M b /M h , where M b =M s + M g and M h is 
the total (dark+baryon) virial mass. 

In Figure [2] V max vs M s (the stellar Tully-Fisher relation, 
sTFR, panel a) and R e vs M s (panel b) at z = are plotted for 
all simulated galaxies (geometric symbols). The black solid 
and dotted lines show linear fits and their estimated la intrin- 
sic scatters to observed normal disk galaxi es in a large mass 
and su rface density range as reported in lAvila-Reese et al.l 
(120081) . Their stellar masses we re corrected by — 01 dex to 
go from the 'diet Sal peter' IMF dBell et~aT1l2003l) implicit in 
their inferences to the IChabrierl (120031) IMF. 

Big and small geometric symbols are used for high- and 
low- resolution runs, respectively. The symbols are open 
when t off = 40 Myr and solid when t Q g = 20 or 10 Myr. Cir- 
cle, pentagon, and square are used for runs where £sf= 0.5, 
0.2, and 0.67, respectively; a triangle is used when the stellar 
mass particle is calculated rather than assigned as a fraction of 
the gas mass cell (see §§2.2). The runs that correspond to the 
same galaxy/halo system have the same color. Symbols tra- 
versed by a cross correspond to high-resolution simulations 
with more grid cells (A* a , B^, and C^). The open inverted 
triangle corresponds to the simulation with the run-away stars 



scheme and without turning off cooling. 

From upper panel of Fig. |2] one clearly sees that the low- 
resolution simulations produce galaxies with too large values 
of Knax for their masses. It is well known that the lack of 
resolution in this kind of simulations introduces an artificial 
dissipation of energy and angular momentum in the gas, pro- 
ducing these too-concentrated galaxies and therefore, peaked 
rotation curves (as well as less SF and energy input to the 
gas, in such a way that galaxy outflows are inefficient, hav- 
ing then the galaxies too high baryonic and stellar fractions, 
see Table 2). However, our high-resolution simulations are 
already in rough agreement with the observationally-inferred 
sTFR, except for the run /. In general, the mid-plane gas 
disk rotation velocity p rofile, V mt (R), is lower than V C (R) in 
the inner regions (e.g.. IValenzuela et al. 1 120071 C10), and in 
some cases, even the maximum of V ro t (R) may be slightly 
lower than V max . Therefore, when comparing with obser- 
vations, Vmax tends to be an upper limit. In Fig. [2] we 
also repr oduce results from othe r recen t high-resolution sim- 
ulatio ns dGovernato et al.l I2007L 1201 Ot iPiontek & Steinmetzl 
1201 lb . where the problem of obtaining a TFR shifted to the 
high-velocity side seems to have been partially overcome (see 
§§4.1 for details on these works). 

Our smallest simulated galaxies fall systematically above 
the extrapolated sTFR of higher masses, suggesting a bend 
towards the low M s (high V max ) side. Such a bend at 
low masses in the sTFR has been discussed recently by 
Ide Rossi et al.l d2010l) . who explain it as a consequence of the 
strong stellar feedback-driven outflows. From the observa- 
tional side, some authors have found that the infrared or stel- 
lar TFRs at low masses become very scattered but apparently 
M a is on average lower for a giv en Vmax than the extrapola- 
tion of the high-mass sTFR (e.g., lMcGaughll2005l : iGeha et alJ 
l2006tlDe Rijcke et alJl2007h . In the upper panel o f Fig. H the 
observ ational data for dwarf galaxies reported in IGeha et aLI 
d2006l) are reproduced (dots). 
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log M. [MJ 

FIG. 2. — Maximum circular velocity (upper panel) and effective radius 
(lower panel) vs. M s for the different galaxies simulated here (figure sym- 
bols, see text for the symbol code and Table 1 for the color/letter code showed 
belo w in lower panel). The sma ll skeletal symbols are for the simulations 
from lPiontek"& Steinmeti EoTTl) and Govemato et al. (2007, 2010). The fits 
and Ict int rinsic scatters in bo t h pane ls to a compiled sample of normal disk 
galaxies by Avila-Reese et al. 12008) are plotted with solid and dotted lines, 
respectively. The observed dwarf galaxies by Geha et al. (2006) are also plot- 
ted (dots). In low er panel, the fit to the blue-cloud disk-dominated galaxies 
from the SDSS by Dut ton et alj I2010bl) is reproduced ( green solid line; the 
dotted lines show the 84th and 16th percentiles of the size distribution). The 
dashed lines are extrapolations to lower masses of the given fits. 

Regarding the R e -M s relation, the lower panel of Fig. [2] 
shows that the low-resolution simulations produce too con- 
centrated galaxies, in agreement with the resolution argu- 
ments already mentioned above. However, even for the high- 
resolution simulations, our most massive modeled galaxies 
have radii smaller than the mean of the observational infer- 
ences for a given M s . In Fig. [2] is also plotted the fit to the 
R e -M s relation inferred recently by iDutton et al.l (1201 Obi) for 
the blue-cloud disk-dominated galaxies from the SDSS (green 
solid line). These inferences show a pronounced bend (flat- 
tening) in the R e -M s relation at low masses ( M s ^ 2 10 10 
M^). The measured radii for dwarf galaxies from lGeha et al.l 
d2006l) are also reproduced in Fig. |2](dots). These galax- 
ies lie c learly below the extrapolation of the the IDutton et al.l 
(1201 0b|) R e -M s relat i on, bu t agree with the extrapolation of 
the lAvila-Reese et alJ (120081) relation. Our lowest-mass simu- 
lations show better agreement with observations. 

It should be not e d that R e both in I Avila-Reese et al.l d2008l) 
and IDutton et al.l (1201 Obi) was calculated assuming R e = 
1.67 Rd, where Rd is the scale radius obtained from fitting an 
exponential law to the disk component of the surface bright- 
ness profile. The above relation is a good approximation 
if most of the stellar mass lies in a disk; however, if the 
bulge/spheroidal component is non-negligible, then the actual 
effective radius R e is expected to be different from 1.67 Rd- 
For several of our high-resolution simulations, we have found 
that R e is smaller than 1.67 Rd by factors of ~ 1.1 — 1.3, not 
enough to explain the differences seen in Fig. |2]for the more 



massive galaxies. In some cases, R e results even larger than 
1.67R d . 

The structure of the stellar component of normal low-mass 
and dwarf galaxies seems to be different from the usual one 
of normal high-mass galaxies, both from the point of view of 
observations and simulations. In particular, the existence of 
an extended spheroidal component could be the rule as seen 
in our lowest-mass simulations (see also C10: lBekkHl2008l 
and therein references on observational works). In general, 
the stellar surface density profiles of our simulated galaxies 
tend to be exponential within ~ 1 — 3 R e . 

Regarding the circular velocity profiles, V C (R), they are 
roughly flat at 1.5 — 3 R e for all the high-resolution runs 
except for the most massive one, Ihi, which shows a nar- 
row peak at 0.2R e . In some cases, V C (R) is even slightly 
increasing at 1 — 2 R e as it is the case of run B* hl , shown at 
z = 0.0 (solid line) and 1.0 (thick dot-dashed line) in Fig- 
ure|3] For our lower-resolution simulations, in some cases the 
velocity profiles are peaked. As expected, V C (R) increases 
with time. The decomposition of V C (R) into the dark mat- 
ter (DM), stellar (s) and gas (g) components is also plotted 
in Fig. [5] for the z = profile. This galaxy is dark-matter- 
dominated for radii >, lR e . The effect of stellar feedback is 
crucial in order to prevent excessive baryon matter concentra- 
tion in the center -and consequently a peaked V c (R) profile- 
as has been discussed in C10 (see alsolGovernato et al]|2007b 
ICeverino & Klvpinll2009l;IStinson et alj|2009h . 

Figure |4] presents the "archeological" SF histories (SFHs) 
for our high-resolution runs, defined as the instantaneous SFR 
as a function of time, and computed for each run using the last 
snapshot recorded. Specifically, in each data dump, in addi- 
tion to the positions and velocities for all stellar particles, the 
code also saves the time at which they formed, their masses 
(initial and present), and metallicities due to SNe of types la 
and II. We add up the (initial) masses of the stellar particles 
formed during a certain time interval, which we take as 0.5 
Gyr, and divide it by this time, to obtain the "instantaneous" 
SFR. Since the amount of stellar mass outside the galaxy is 
typically a small fraction, the SFH computed inside i? v j r is 
mostly that of the central galaxy. Unlike C10, where a time in- 
terval of 0. 1 Gyr was taken to compute the SFHs, here we use 
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FIG. 4. — "Archaelogical" SFHs from the high-resolution runs. Top-panels: 
Ihl> Fhli Ehl mns (left panel) and D^i, A^i rans (right panel). Lower- 
panels: A comparison of the B ran with and without "*"; that is, run with a 
low and high number of cells inside -R v i r (left panel). Right-panel: Same as 
in left panel but for ran C. With the exception of ran F^i, models show a 
much smaller SFR in the last 5-6 or more Gyr with respect to the previous, 
active, early epoch. 

0.5 Gyr in order to produce smoother SFHs. To estimate how 
much the SFR fluctuates, for example, in the last 1 Gyr, we 
measured the SFR in bins of 0. 1 Gyr width and computed the 
dispersion of data for each run. We find dispersions around 
the means of about 30-40%. 

In the top panels, the "archeological" SFHs of runs 
I hl ,F hl ,E hl (left panel), and D hl ,A hl (right panel), are 
plotted, using the same color coding as in Figure 1 . In the 
lower-left panel, we compare the SFHs of run B simulated 
with a high number of DM matter particles but with a low 
(dotted blue lines) and high (solid black lines) number of 
cells. The same comparison but for run C is shown in the 
lower-right panel. As can be seen, there is not a significant 
difference in their SFHs between models with "*" and those 
without it. The SFHs are diverse but in most cases they are 
characterized by an active phase that starts after the first 1 -2 
Gyr and lasts 2-6 Gyr. This phase is followed by a "quies- 
cent" stage with a low value of the SFR, with most of the SF 
concentrated in the central part of galaxies. In most cases, 
the SFHs inherit some features of the halo MAHs. The "ar- 
chaeological" SFHs already forecast that simulated low-mass 
galaxies will have low SSFRs at late epochs in conflict with 
observational inferences. 

Finally, it should be said that in almost all of our simulations 
a multi-phase ISM develops in the disks (see for more details 
C10). This results from the efficiency by which thermal en- 
ergy is injected into the ISM and the ability of the medium to 
self-regulate its SF. 

3.2. Specific star formation rate as a function of mass 

The main result of this paper is shown in Figure where 
the directly measured galaxy SSFR of the simulations at four 
epochs (z = 0.00, 0.33, 1.00 and 1.50) are plotted vs the cor- 
responding current stellar masses (geometric symbols). For 
the symbol code and color/letter code, see previous subsec- 
tion. Unlike the "archeological" SFR shown in Fig. [4] here 



the SFR at a given redshift is computed by summing up the 
mass of all stellar particles with ages smaller than 0. 1 Gyr lo- 
cated in a cylinder with 1.0 /i _1 kpc (proper) height and 20.0 
/i _1 kpc of radius, centered on the gaseous disk. The SFR 
is then this mass divided by 0. 1 Gyr. This time bin is wide 
enough so as to avoid short fluctuations in the SFR values. 

The general trend of the measured SSFR in the simulations 
is to decrease with time, though, in general the SFR (and 
therefore SSFR) histories are quite episodic (see also Fig. 0J. 
It is remarkable that at low redshifts (z = and 0.33), all sim- 
ulated galaxies have SSFRs below the SSFR a galaxy would 
have if it had formed all stars at a constant SFR 7 (horizontal 
dashed line). This means that the M s of these galaxies was as- 
sembled relatively early, the current SFRs being smaller than 
their past average SFRs (the galaxies entered into a quiescent 
SF phase). For higher redshifts, z = 1.00 and 1.50, the SS- 
FRs of several of the high-resolution simulations tend to be 
already around the SFR=const. line, which means that galax- 
ies are in their active phase of SF and M s assembly. 

Notice that some runs show negligible SSFRs at some 
epochs. For example, the SSFR of run At 1 is zero at z = 0.0 
(the halo MAH of this run is the earliest one among all simula- 
tions, see Fig. [T] which seems to imply that its SFR decreased 
strongly very early, see Fig. |4). Because of the episodic na- 
ture of the SFR it may well be that, even though we use a 
relatively wide time bin to compute the instantaneous SFR 
(0. 1 Gyr), there are some epochs for which we see a galaxy 
almost not forming stars. Run Cn is special in the sense that 
it has zero SFR at the time-steps corresponding to z = and 
z = 0.33, something that is not seen for the higher-resolution 
counterpart simulations. The resolution in this case seems to 
affect the efficiency of SF because the galaxy actually is grow- 
ing plenty of gas (f g = 0.85 and 0.70 at z = and 0.33, 
respectively; note that the halo MAH of this model shows 
the most active growth at late epochs). This is probably be- 
cause instabilities in the disk have not been resolved. Indeed, 
in the corresponding high-resolution runs, specially the 
one, a relative high SSFR is measured at late epochs. How- 
ever, there are cases of high-resolution simulations, for exam- 
ple run Dhi, where the gas fraction is high and even increases 
at late epochs, but the SSFR keeps low and decreasing with 
time. We suspect that in these cases some physical process is 
damping the disk gas instabilities, and therefore inhibiting ac- 
tive SF; the presence of a dynamically hot spheroid -which is 
actually formed in several of our lo w-mass simulations- has 
been shown to work in this direction dMartig et al.ll2009l) . 

From Fig. |5]one sees that the resolution affects moderately 
the SSFRs. Our low-resolution simulations tend to have lower 
values of SSFR at z = 1.50 and 1.00 than the high-resolution 
ones, but the opposite applies for lower z. At this point, it 
should be said that it is difficult to establish any systematicity 
due to the episodic behavior of the SFR in the simulations. At 
2 = 0, low-resolution simulations seem to produce slightly 
higher SSFR values than high-resolution simulations. Note 
also that the former have systematically higher values of f g 
than the latter (see Table 2). The effect of increasing the num- 
ber of cells (filled circles traversed by a cross) is not so signif- 
icant. The global effect of the strength of the stellar feedback 
(regulated in our simulations by esf and t g) apparently is to 

7 If SFR=const. in time, then SFR/M S = 1/[(1 -R)(t H (z) - 1 Gyr)], where 
R = 0.4 is the average gas return factor due to stellar mass loss, tfj is the 
cosmic time, and 1 Gyr is subtracted in order to take into account the onset 
of galaxy formation. 
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FIG. 5. — SSFR vs M a measured at four different redshifts for all of our simulations (geometrical symbols). A given galaxy is represented with the same 
color (associated to the letter that identifies the simulation according to Table 1), and the different symbols are to distinguish runs of the same galaxy but with 
different sub-grid parameters (see §§.. for the symbol code). The horizontal dotted lines in each panel correspond to the SSFR of a galaxy that had a constant 
SFR in the past (see text); galaxies above (below) this line are expected to have been les s (more) active in forming sta rs on average in the past t han in the current 
epoch. In the z = panel results from si mulations by Gove rnato et al. (2007, 2010), Piontek & Steinmetz 12011), and Sawala et al. |20TT]) are also plotted; 
the solid line is the linear fit carried out by Salim et al. (2007) to the sub-sample of star-forming galaxies from a volume- corrected sample of ~ 50, 000 SDSS 
galaxies (the intrinsic scatter is shown with dotted lines), while the long-dashed line is the Schechter fit to the entire sample, i.e. including composite galaxies 
(SF/AGN), galaxies dominated by AGN emission, and those with no Ho detection. The skeletal symbols, fit lines, and horizontal segments in the high-redshift 
panels correspond to different observational estimates down to the completeness limits reported in the works indicated inside the panels (see Appendix for a 
description of each one of them). At all the epochs reported here, the observationally estimated SSFRs of low mass-galaxies are on average significantly higher 
than those of simulations. While observational estimates at each one of the redshifts are typically above the SFR=const. lines for low masses, the simulations are 
in most cases below them, specially for z < 1. 



keep the SF more active at later epochs as seen for runs E^i 
(t s = 10 Myr, open green circle) and £7,2 (i ff= 40 Myr, 
filled green circle). However, the differences are small, as in 
the case of variations in esF from 0.5 to 0.2 (runs Eya, filled 
green circle, and Em, filled green pentagon, respectively). In 
general, variations of the sub-grid parameters around reason- 
able values seem not to be a source of significant and system- 
atical changes in the SSFR of the simulated galaxies. 

3.2.1. Comparison with observational inferences 

In each panel of Fig. |5j different observational inferences 
compiled here are plotted for comparison. Unless otherwise 
stated, we plot only the data that obey the completeness limit 
given by the different authors. The inferences of M s (and 
SFR) from the observed luminosities or spectral energy dis- 
tribution of galaxies are based on results from SPS (stellar 
population synthesis) models. For all cases reproduced in Fig. 
[5] a constant, universal, stellar initial mass function (IMF) has 
been assumed in these models, although different authors may 
have used different IMFs. W e homogenize all observational 
results to a IChabrierl (120031) IMF (see Appendix for the dif- 
ferent corrections applied here). A description of each one 
of the observational sources compiled from the literature is 
presented in the Appendix. 

From a visual inspection of Fig. we conclude that all 
simulated low-mass galaxies at different resolutions and sub- 
grid parameters lie significantly below the averages of a large 
body of SSFRs observational determinations as a function of 
M s out to z ~ 1.5, especially for the blue star-forming ones. 
In the z = panel, a compilation of recent numerical re- 
sults from other authors are also plotted (skeletal symbols; 
see §§4.1 for more details). 

By taking into account only our high-resolution simulations 
with esF = 0.5, we estimate that the SSFRs in the mass range 



M s ~ 10 s - 10 10 M are - 5 - 10 times lower than the 
means inferred from observations at z w 0. These differ- 
ences are roughly the same or slightly smaller at z w 0.33. 
At higher redshifts (z ps 1.0 — 1.5), the comparison becomes 
difficult by the incompleteness in the observational samples. 
However, both for the small mass range where our simulations 
coincide with observations above the completeness limit and 
the extrapolations to lower masses of the inferred SSFR-M S 
relations, the differences persist, though they are less dramatic 
than at lower z. 

It is notable the systematical shift of the SSFR-Af, rela- 
tion of observed galaxies as z decreases: from a rough av- 
erage of all relations, the typical stellar mass, M quenc h, that 
crosses below the line corresponding to the constant SFR case 
at each z (horizontal short-dashed lines in Fig. |5) decreases 
with time. In other words, the lower the z, the smaller on aver- 
age are the galaxies that start quenching their SF (downsizing 
in SSFR). Interesting enough, the downsizing of M quonc h is 
in rough quantitative agreement with the downsizing of the 
transition mass from active to passive galaxies found from 
the semi-empirical M s growth tracks in FA10: logAftran(z) 
= 10.30 + 0.55z. For the simulated galaxies, there is no evi- 
dence of such a phenomenon of downsizing. 

3.3. Stellar mass fraction as a function of mass 

In Fig. [6] the stellar mass fraction, F s = M s jMh, of our 
simulated galaxies vs Mh at z = 0.00, 0.33, 1.00 and 1.50 are 
plotted. Results from numerical simulations of other authors 
shown in Fig. |5]are also plotted (skeletal symbols). 

The bell-shaped curve in each panel corresponds to the 
continuous analytical approximations given in FA 10 to the 
semi-empirical determinations of the M s -Mh relations at dif- 
ferent redshifts (0 ^ z < 4) performed by iBehroozi et al.l 
d2Q 1 Ot) . These inferences are based on the abundance match- 
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FIG. 6. — Stellar mass fraction vs of the simulated galaxies at z = 0, 0.33, 1, and 1.5 (geometrical symbols). The symbol and color/letter codes are as in 
Fig|5] The results from other authors (indicated insi de the panel) are also reproduced for 2 = 0. The solid curves in each panel are semi-empirical inferences 
as reported in FA10 based on[Behroozi et al. (2010). The error bars show the lcr uncertainty of these inferences. The dotted lines are extrapolations to lower 
masses, given the lower limits in M s of the empirical galaxy M s functions used in Behroozi et al. For the dashed line in the 2 = panel see the text. 



ing formalism, where the cumulative observed galaxy stellar 
mass function at a given epoch is matched to the cumula- 
tive halo mass function in order to find the M/j correspond- 
ing to a given M s , under the assumption of a one to one 
galaxy-halo correspondence. Several authors have found sim- 
ilar results out to z ~ 1, even when different observational 
data sets and formalisms have been used (see for recent re- 
sults e.g., iMoster et all 1201 (A iGuo etal.ll20Tot IWang & Jingl 
120101) . For the local universe, direct techniques, as galaxy- 
galaxy weak lensing and satellite kinematics, have also been 
applied to infer the M s -Mh relation with results that are in 
reasonable agreemen t among them and with other techniques 
(see for comparisons IBehroozi et al.ll2.Q10t iMore et all 1201 Pi 
I Rodrfguez-Puebla et al. 11201 ll and more references therein) 
in the mass ranges where they can be compared (the direct 
techniques are reliable for the time being only for intermedi- 
ate or massive galaxies). The error bar in Fig. |6]indicates the 
approximate ±lcr statistical and sample variance uncertain- 
ties in the inferences. The main contributor to this uncertainty 
com es from assumptions i n converting galaxy luminosity into 
M s (IBehroozi et al.ll2010t) . 

According to Fig. [6] the smaller is the mass, the lower is F s 
on average, both for our simulations and the semi-empirical 
inferences. However, the simulation results have F s values 
systematically higher than those inferred from observations. 
At z ps 0, the differences amount to factors around 5-10 
and they apparently increase at higher redshifts; at z = 1.5, 
the simulations have values of F s for a given M s around 
30 times higher than those inferred semi-empirically at the 
masses where the comparison can be done; the same factor 
applies for the lowest masses if one extrapolates the semi- 
empirical determination to these masses. For simulated galax- 
ies, F s changes very little in the redshift range reported here 
while halo masses decrease as z increases. Therefore one ex- 
pects the F s -Mh relation to shift on average to the low-mass 
side as z gets higher. This behavior is contrary to the semi- 
empirical inferences. 

In the semi-empirical inferences, the mass function of pure 
dark matter halos is used. It could be that the masses of the 
halos (dark+baryonic matter) end up smaller when baryonic 
processes are included; for example, as a result of mass loss 



out the virial radius in low-mass halos. The maximum frac- 
tion of ejected baryons at z = in our simulations is 82% for 
runs A* tl and i?/a 8 - This would imply that Mh in run 
(or B h i) is ps (1 - OMF^u) -1 = 1.14 times smaller than 
in a purely dark matter simulation. The dashed line in Fig. [6] 
z = panel, shows the expected shift in the M s /Mh — Mh 
relation by this maximum factor. On the other hand, the in- 
ferred MJMh ratio at low masses depends mainly on what 
is the faint-end slope of the galaxy M s function. If observa- 
tions show that there is a steepening of this slope, then the 
M s IMh ratio may flatten at low masses, being in better agree- 
ment with our simulation results. Some recent observational 
studies suggest that the higher the z (for z £ 1), the steeper 
is the faint-end slope dKaiisawa et alj 12009; Mo rtlock et alj 
1201 ll) . For lower z, t here is also eviden ce of a somewhat steep 
slope (~ —1.7. e.g. iDrorv et alj|2009l) . However, even that, 
the MJMh-Mh relation at low-masses seems t o be steeper 
and lo wer than in our simulations (see Fig. 9 in IDrorv et aD 
l2009h . 

In Fig. [7] we have merged into one the four panels (epochs) 
of Fig. [6] showing the evolution (connected with dotted lines) 
of only some of our high-resolution simulations. The rela- 
tively large jumps seen in some runs are mainly due to halo 
major mergers. For example, run Fhi shows a strong decreas- 
ing of F s =M s IMh at z ~ 1 because it suffered major mergers 
around this epoch (see Fig. [!)■ For simulations where the halo 
MAH tends to be smooth, F s changes more smoothly with z 
(e.g., runs D hl , E hl and I hl ). 

We also plot in Fig. [7] the evolution of two semi-numerical 
models presented in FA 10. These models correspond to the 
case of energy-driven SN outflows and no re-accretion, with 
an outflow efficiency of 65%; only with such a high efficiency 
it is possible to lower the stellar mass fractions to roughly 
agree with the F s -Mh determinations at z = 0. 

For both cases, simulations and semi-numerical models, F s 
changes relatively li ttle with z. As mentioned i n the Introduc- 
tion, FA10 (see also lConrov & Wechslerll2009l) connected the 
F s -Mh relations by using average MAHs and obtained the 

8 For other mass runs with Mh < 3 X 10 9 Mq, this fraction vary around 
60-80%, while for the more massive runs the fractions are 10 — 40%. 
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FIG. 7. — Stellar mass fraction vs at z = 0, 0.33, 1.0 and 1.5 for the 
high-resolution runs indicated inside the panel, and for two semi-numerical 
models presented in FA 10 (open black triangles; the dotted lines extend up 
to 2 = 4). The solid black curves are the same semi-empirical determina- 
tions plotted in Fig. [6] at the four epochs mentioned above. The dashed red 
curves correspond to three semi-empirical individual average evolutionary 
tracks (GHETs) as inferred in FA 10. 

corresponding individual (average) F s (or M s ) hybrid evolu- 
tionary tracks (GHETs, dashed red curves in Fig. |7). These 
tracks in the logM s (z)(or logF s )-logA J f/ i (z) plane have very 
steep slopes for low-mass galaxies, dlogAI s (z)/dlogMh(z)~ 
4.5 (z < 1), i.e. for these galaxies, M s (F s ) grows very fast 
since z ~ 1. However, for the semi-numeric ally modeled 
galaxies, d\ogM s (z)/d\ogM h (z)~ 1.4, i.e. M s (F s ) grows 
slowly. Our numerical simulations, albeit noisily, confirm 
such a behaviour, which strongly disagrees with the semi- 
empirical tracks. 

4. DISCUSSION 

We would like to stress that (i) the potential issues of low- 
mass ACDM-based simulated galaxies showed in Figs. [5]and 
[6] seem to be shared by most previous high-resolution numeri- 
cal simulation works, and (ii) that the empirical inferences we 
have compiled in this paper are yet subject to large uncertain- 
ties and probable systematical errors. Following, we discuss 
both items. 

4. 1 . Previous numerical works and numerical issues 

In Figs. [2] [5] and [6] we have reproduced the z = proper- 
ties of sub-M* galaxies presented in recent numerical works, 
where very high-resolution re-simulations of individual low- 
mass galaxies were carried out and the quantities we are in- 
terested in were reported; 

The simulations by iGovernato et al.l (120071 1201 Ol) (aster- 
isks in Figs. 12] [5] and|6]l were performed with GASOLINE, 
a N-body + Smoothed Particle Hydrodynamics (SPH) code. 
The SF prescription ensures that the SFR density is a func- 
tion of the gas density according to the observed slope of the 
Kennicutt-Schmidt law, and a SF efficiency parameter sets the 
normalization of this relation. In order to allow for SN-driven 
expanding hot bubbles in the non-resolved regions, cooling is 
turned off in the gas particles receiving SN energy until the 
end of the snowplow phase according to the Sedov-Taylor so- 
lution. For the force resolutio ns of ~ 100 — 200 pc a chieved 
in their low-mass simulations, IGovernato et al.l (1201 Ol) used a 



high value for the SF density t hreshold, usf = 100 cm " 3 (for 
the high-mass simulations in IGovernato et all 120071 ngp = 
0.1 cm~ 3 was used), in order to have enhanced gas outflows 
that remove low angular momentum gas from the central re- 
gions of the galaxy; according to the authors, this helps to 
create realistic disks. 

The simulations bv lPiontek & Steinmetzl d201 lb (crosses in 
Figs. |2 [5] and [6]l were performed with GADGET-2, a Tree- 
PM N-body + SPH code. The SFR was implemented accord- 
ing to a Schmidt law, /?* = c*p g /tsF, where c* is a SF ef- 
ficiency parameter, and a stochastic approach for assigning 
stellar particles to the gas particles masses was used. Cooling 
is turned off in the gas particles receiving SN energy by a fixed 
time, 20 Myr. We reproduce here their "standard model" high- 
resolution simulations (gravitational softening parameters of 
100— 140 pc). With few exceptions, they find little or no angu- 
lar momentum deficiency in their galaxies, though prominent 
bulges persist in most cases and galaxies lie slightly shifted in 
the stellar TFR. 

The simulations by ISawala et al.l (1201 ll) (skeletal triangles 
in Figs. [5] and |6]l correspond to six representative low-mass 
halos (M h (z = 0)~ 10 10 M Q ) extracted from the Millenium- 
II Simulation. These halos were resimulated at high reso- 
lution with gas (dark and stellar mass particles of 8 x 10 4 
and 1.8 x 10 4 M©, respectively) using the Tree-PM code 
GADGET-3, which includes metal-dependent cooling, SF 
(Schimdt law and stochastic approach), chemical enrichment, 
and energy injection from type II a nd la SNe implemented in a 
multiphase gas model developed in lScannapieco et aTI(l2006l) . 
SN energy is shared equally between the hot and cold phases. 
Cold particles which receive SN feedback accumulate energy 
until their thermodynamic properties raise to the typical prop- 
erties of the local hot phase. The gravitational softening scale 
within the collapsed halo is 155 pc. The final simulated ob- 
jects have structures and stellar populations consistent with 
observed dwarf galaxies 

On one hand, as seen from Figs. [2] |5] and |6] in most cases 
the properties of simulated (central) low-mass galaxies for a 
given M s , in particular SSFR and F s , agree among different 
authors including our results, in spite of the different codes 
used, simulations settings, and sub-grid physics, as well as 
the diversity in halo MAHs. On the other hand, all these simu- 
lated galaxies tend to have realistic structural, dynamical, and 
gaseous properties, excepting the too high SSFR and too low 
F s values at a given M s as compared with observational in- 
ferences out to z ~ 1. 

As said in the Introduction, by increasing the resolution and 
including adequate SF/feedback prescriptions, several of the 
difficulties found in previous numerical simulations of disk 
galaxies have been overcome, specially the one related to the 
angular momentum catastrophe. According to simulation re- 
sults, a key physical ingredient in the evolution of disk galax- 
ies is the SF-driven feedback (see discussions in C10 and in 
the above mentioned references): it self -regulates SF promot- 
ing a multi-phase ISM; it removes low angular momentum 
gas from the galaxy avoiding this way the formation of a too 
compact stellar component with very old stellar populations; 
it drives galaxy outflows that lower the baryonic and stellar 
mass fractions in the simulated galaxies. Most authors agree 
that without a highly efficient SF-feedback, realistic galaxies 
are not reproduced. In our simulations, such an efficiency 
is attained by artificially turning off the cooling by a time 
toff (only) in the cells where stars form (§§2.1). The times 
t a g used here (10-40 Myr) correspond roughly to the time a 
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pressure-driven super-shell takes to reach £ 100 — 150 pc. 
This latter scale corresponds roughly to our spatial resolution 
limits. Tests show that for t g varying in the 10-40 Myr range, 
our results do not change significantly; compare, for example, 
runs Ehi and Eh2- 

Regarding SF, the density and efficiency parameters used 
in the SF schemes, sh owed to be dependent on the resolution 
of the simulation (e.g JSaitoh et alJl2008l;lBrooks et afll201 lb 
and interconnected with the feedback and gas infall processes. 
From physical considerations (see C10), for the spatial resolu- 
tions achieved here, the value for nsF we ad opted is 6 cm" 3 . 
If tisf is increased to 50 — 100 cm -3 as in iGovernato et al.l 
(1201 Oh . then our low-mass galaxies become too concentrated 
and having too early SF (see also C10). 

The lack of resolution at the scales where SF and the SF- 
driven momentum and energy injection to the ISM happen, 
remains as a shortcoming of numerical simulations, in spite 
of the mentioned improvements in the implementation of 
sub-grid schemes. The sub-grid schemes are typically "op- 
timized" to reproduce present-day galaxy properties and for 
following the highest density regions of the simulation, where 
galaxies form. Less attention has been given to the physics of 
the intra-halo and intergalactic medium. The processes re- 
lated to this medium could play a relevant role in galaxy for- 
mation and evolution, in particular at the early epochs and for 
low-mass halos. 

It is not clear whether the main physics and evolutionary 
process of disk galaxy evolution are already captured by our 
and other simulations and the challenge remains only in im- 
proving reso lution and tuning better the combination of small 
effects (e.g. JPiontek & Steinmetzl201 lb or some key physical 
ingredients and evolutive processes are yet missed. The an- 
swer to this question is beyond the scope of this paper. How- 
ever, the results and discussion presented here point out to a 
potential serious shortcoming of simulations and models of 
low-mass central galaxies in the context of the ACDM sce- 
nario. The solution to such a problem requires that the stellar 
mass assembly of galaxies in distinct low-mass halos be sig- 
nificantly delayed to late epochs, and the smaller the galaxy, 
the more should be such a delay (downsizing). 

4.2. Observational caveats 

In the last years we have seen an explosion of works aimed 
to determine the stellar mass and SFR of galaxies at differ- 
ent redshifts (see Introduction and Appendix for references). 
Yet these inferences should be taken with caution, though the 
problem stated in this paper seems to be robust. Following, 
we discuss the main concerns about these inferences. 

4.2.1. Uncertainties in the inferences of M s and SFR 

The SPS models used to fit observational data and infer 
hence M s of galaxies are flawed mainly by uncertainties 
in stellar evolution (for example in the thermally-pulsating 
asymptotic giant branch, TP-AGB, and horizontal branch 
phases) and by the poor knowledge of the initial mass func- 
tion, IMF, as well as due to degeneracies like the one be- 
tween age and metallicity (see for recent extensive discussions 
Maraston et al. 2006; Bruzual 2007; Tonini et al. 2009; Con- 
roy, Gunn & White 2009; Santini et al. 2009; Salimbeni et 
al. 2009). For example, Conroy et al. (2009) estimated that 
including uncertainties in stellar evolution, M s at z ~ carry 
errors of ~ 0.3 dex at 95% confidence level with little depen- 
dence on luminosity or color. 



Maraston et al. (2006) showed that the stellar masses of 
galaxies with dominating stellar populations of ~ 1 Gyr age 
could be on average 60% lower if their assumptions for con- 
vective oversh ooting during the TP-A GB phases are used. 
More recently, ISalimbeni et al.l (120091) have determined the 
stellar masses of galaxies from a GOODS-MUSIC sample 
at redshifts 0.4 < z < 4.0 by using the Bruzual & Char- 
lot (2003, BC03) and Chariot & Br uzual 2007 (CB0 7, see 
Bruzual 2 007) SPS model s and the iMarastonl (120051) mod- 
els (M05). ISalimbeni et all d2009l) have found that the masses 
obtained with MR05 and CB07 (which take into account the 
TP-AGB phase) are on average lower than those obtained 
with BC03. In the redshift bin 0.4 < z < 0.8 the ratio they 
have found between BC03 and BC07 masses is ~ 0.3 dex for 
M s < 10 9 - 5 Mq and ~ 0.05 dex for M s > 10 11 M Q , while 
between BC03 and M05 the ratio is ~ 0.2 dex for all masses. 

All the SSFR-M S relations compiled here and showed in 
Fig. [5] were inferred using the old BC03 SPS models. If we 
assume a conservative correction due to the TP-AGB phase 
of 0.2 dex in M s for all the masses and a slope of —0.5 
in the logSSFR-logA/ s relation (in most cases this slope is 
shallower), then such a relation shifts on average by 0. 1 dex 
(« 25%) towards the high-SSFR side. Therefore, if any, the 
differences between observed and predicted SSFR-M S rela- 
tions could be even larger than showed in Fig. [5] 

Regarding the determination of SFR, it coul d significantly 
depen d on the used indicator. For example, in iGilbank et al.l 
d201 ll) is showed that the SFRs inferred for a sample of galax- 
ies at z ~ 1 by using as a indicator the [Oi l] empirically 
correc ted for extinction as a function of M s (IGilbank et al.l 
1201 Ol) are on average 2.2 times lower than those inferred by 
using [Oil] + 24/im, where the [Oil] has not been corrected 
for extinction (since the light from SF being reprocessed by 
dust should now be measured by the 24/im flux). This ex- 
plains the difference see n in Fig. [5] thi r d panel , in between 
IGilbank et all ( 120111) and lNoeske et al. I d2007bl) (and proba- 
bly the other authors). 

Different indicators have different sensitivities to unob- 
scured/obscured SF. Because observations are required at dif- 
ferent z and for different mass and galaxy type ranges, dif- 
ferent indicators in the same sample are often used. This un- 
avoidably introduces systematics that difficult the interpreta- 
tion of correlations and evolutive trends. It should be said that 
in several of the works compiled here, the combination of two 
or more indicators that trace both obscured/unobscured SFR 
were used with calibrations given by SPS models. These cali- 
brations seem to be robust to variations in the stellar evolution 
and extinction assumptions, but they vary among different as- 
sumption on IMF by w 0.1 — 0.4 dex. The latter system- 
atic variations are partially compensated when calculating the 
SSFR because both the SFR calibrators and the estimated M s 
change in the same direction with varying the IMF. Besides, 
the effects of extinction are actually minimal at low masses. 

In the compilation of observational determinations showed 
in Fig. [5J a large diversity of SFR indicators, extinction cor- 
rections, and galaxy samples were used. Although there are 
differences in the results among these works, our conclusions 
do not change when using one or another particular work. 
This means that among all the current observational determi- 
nations, there is a rough convergence regarding the SSFR- M s 
relation and its evolution. 

4.2.2. Selection effects 
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Other significant sources of uncertainty and possible sys- 
tematics in the inferred SSFR- Ms relations at different z are 
(e.g., Daddi et al. 2007; Chen et al. 2009; Stringer et al. 
2010): selection effects due to sample incompleteness in a 
given wavelength; environmental effects; limit detections of 
the indicators of SFR due to flux-limits or low emission-line 
signal-to-noise ratios; in addition, obscured AGN emission 
could contaminate the infrared flux and some of the optical 
lines used to estimate SFR, though at low masses AGNs it is 
not frequent. In most of the works presented here, the authors 
sought to refine their samples to minimize the above selection 
effects and the contamination by AGNs. 

The main concern among the mentioned above issues is that 
selection and environmental effects could bias the observed 
SSFR-M S relations. For example, the high SSFRs of low- 
mass galaxies could be due to transient star bursts if the SF 
regime of these galaxies is dominated by episodic processes. 
If among the low-mass galaxies those with low SFRs are 
missed due to detection limits, then the SSFRs of low-mass 
galaxies will be biased on average toward higher SSFRs, a 
bias that increases for samples at higher redshifts. Neverthe- 
less, in most of the observational studies reporting the SSFR- 
M s relation, the authors discuss that while this is possible at 
some level, it would hardly change the tight S SFR- Ms r ela- 
tions observed since z ~ 1 (see e.g jNoeske et al.ll2007al) . 

On the other hand, the analysis of very local surveys helps 
to disentangle whether episodic star bursting events dominate 
or not the SF history of low-mass galaxies 9 . From a study 
of the SF activity of galaxies within the 1 1 Mpc Local Vol- 
ume, Lee et al. (2007; see also Bothwell, Kennicutt & Lee 
2009) have found that intermediate-luminosity disc galaxies 
(-19 <Af B < - 15 or 50 < V max / km s" 1 < 120) show 
relatively low scatter in their SF activity, implying factors not 
larger than 2-3 fluctuations in their SFRs; above Vmax ~ 120 
km/s the sequence turns off toward lower levels of SSFRs 
and larger bulge-to-disc ratios. These results are for nearby 
galaxies, where selection effects are minimal, and imply that 
the SSFRs of disc galaxies with M s > 5 10 s M Q follow a rel- 
atively tight sequence, without strong fluctuations. For galax- 
ies smaller than V max ~ 50 km/s (dwarfs) the situation seems 
different. The results by Lee et al. (2007) show that a signif- 
icant fraction of such galaxies are undergoing strong episodic 
SF fluctuations due to the large scatter in their SSFRs. 

Another observational study of nearby galaxies by James 
et al. (2008), also concluded that there is little evidence in 
their sample of predominantly isolated field galaxies of sig- 
nificant SF through brief but intense star-burst phases. There- 
fore, it seems that the tight sequence found for normal star- 
forming galaxies in the SSFR-M S plane in large surveys as 
SDSS (Brinchman et al. 2004; Salim et al. 2007; Schimi- 
novich et al. 2007) is intrinsic and due to a high degree of 
temporal self-regulated SF within indi vidual galaxies. This 
sequence (called the 'main sequence' in lNoeske et al.ll2007al) 
seems to persist back to z ~ 1 as discussed above. 

Finally, possible selection biases in the SSFR-M S relation 
at low masses due to environment seem not to be a concern. 
For local SDSS and high z (up to z ~ 1) zCOSMOS galaxy 
samples it was found that the relationship between SSFR and 
M s is nearly the same in the hig hest and lowest d ensity quar- 
tiles of star-forming galaxies dPeng et al.l l2010h . What is 
strongly dependent on environment are the fractions of star- 

9 Have in mind that blue, late-type galaxies are the dominant population in 
the mass range 10 9 < M s /M < 5 10 10 M Q (e.g. lYang et al.120091) . 



forming and quenched galaxies: the higher the environmen- 
tal density, the higher the fraction of quenched (red) galax- 
ies, even at lower mass. On the average -where most galax- 
ies live- and low density environments, the fraction of red 
galaxies becomes smaller and smaller as the ma sses are lower 
( Peng et alJ2010h . Similar results were found by iMcGee et al.l 
d2010l) . who determined stellar masses and SFRs for large 
samples of field and group galaxies at z ~ by using SDDS 
and at z ~ 0.4 by using the Galaxy Environment Evolution 
Collaboration survey. The SSFR- Ms relation of star-forming 
galaxies is consistent within the errors in the field and group 
environment at fixed z, but the fraction of passive (red) galax- 
ies is larger in groups than in the field at almost all masses. 

4.2.3. Concluding remarks 

We conclude that, in spite of significant uncertainties and 
possible systematics due to selection effects, current observa- 
tional determinations of the SSFR-M S relation for field galax- 
ies at redshifts up to z ~ 1 have achieved a rough consistency 
among them. The global empirical picture of the M s assem- 
bly that they suggest, for sub-M* central galaxies, begins to 
be established; however, several important details such as, the 
separation by galaxy types and environment, are still highly 
uncertain. According to those works, where the samples 
were separated into blue/star- forming, red/quenched and/ or 
AGN- dominatd galaxies (e.g., ISalim et all 120071: iBell et al. I 
l2007t iNoeske et ai1l2007al: iKarimet al.ll2011l) . the low-mass 
side of the SSFR— M s relation is actually dominated by 
blue/star-forming galaxies, which correspond typically to 
disk-dominated galaxies. 

Big efforts should be done in the next years to determine 
with better accuracy and completeness the M s and SFR of 
galaxies down to small masses, up to high redshifts, and sep- 
arated by galaxy type and environment. This implies not only 
observational efforts but also theoretical ones: SPS models, 
used to infer the physical quantities from the flux and spec- 
troscopy measurements, need to be improved. The informa- 
tion provided by the SSFR as a function of M s , at different 
epochs, offers important clues for understanding how galaxies 
assembled their masses and helps to constrain any theoretical 
approach to galaxy formation and evolution. 

5. CONCLUSIONS 

The N-body + Hydrodynamic ART code has been used to 
simulate sub-M* central galaxies in the halo mass range of 
M h (z = 0)— 2 - 40 x 10 10 M Q , varying the resolution 
and several of the sub-grid parameters. Most of the obtained 
galaxies in our highest resolution simulations have dynami- 
cal and structural properties in reasonable agreement with lo- 
cal sub-M* field galaxies, which typically are disk, late-type 
galaxies (Figs. |2]and[3] and TableO. 

Although we have simulated only 10 different galaxies, all 
immersed in different distinct halos, our results allows us to 
answer, at least preliminary, the main questions stated in the 
Introduction and hence conclude that: 

• The SSFRs at z ~ and ~ 0.3 of simulated galaxies in 
the mass range M s ~ 10 8 - 10 10 M© are 5 - 10 times 
lower than the mean values inferred from several obser- 
vational samples of field galaxies (Fig. 0. There is not 
a high-resolution simulation with a SSFR value high 
enough to lie above the lo wer la scatter of th e obser- 
vational determinations by ISalim et ail (120071) of star- 
forming galaxies at z <~ (most observed low-mass 
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central galaxies are actually star-forming). At higher 
redshifts, z ~ 1.0 — 1.5, most simulated galaxies have 
masses below the completeness limits of current ob- 
servational inferences at those z. However, both the 
measurements of the SSFRs of those galaxies which 
are below these limits and the extrapolations to lower 
masses of the SSFR-A/ S relations of the complete sam- 
ples, are higher than the SSFRs of our simulated galax- 
ies, though the differences are apparently smaller than 
at z ~ 0. Several of the simulations at z ~ 1.0 — 1.5 
have already SSFRs around the value corresponding to 
a galaxy forming stars with a rate equal to its past aver- 
age. On the contrary, at z ^ 0.3, most simulations have 
SFR values well below this case, in clear disagreement 
with most current observational determinations. 

• The evolution of the observationally determined SSFR- 
M s relations shows a clear trend of downsizing for the 
typical stellar mass, Mq Uenc h, of galaxies that transit 
from active to passive (those that cross below the line 
of constant SFR at a given z in Fig. |5J. There is not any 
evidence of such a downsizing trend in the simulations. 

• The stellar mass fractions {F s = M 8 /Mh) of simulated 
galaxies are 5 — 10 times larger at z ~ than current de- 
terminations (semi-empirical and direct) of these frac- 
tions as a function of M s (Fig. [6). At higher redshifts, 
the differences, at a given mass, increase even more; at 
z = 1.5, the F s of simulated galaxies are around 30 
times larger than the semi-empirical inferences or their 
extrapolations to lower masses. Put in another way, 
while the F s -Mh relation of simulated galaxies would 
tend to remain constant or slightly shift to higher F s 
values as z increases, the current semi-empirical infer- 
ences show that the low-mass side of this relation shifts 
significantly to lower F s values (see also Fig. |7). 

Therefore, our numerical simulations confirm the issues 
previously found with semi-numerical models of disk galaxy 
evolution (and eventually SAMs) in the context of the ACDM 



scenario: low-mass disk-like galaxy models seem to assemble 
their stellar masses on average much earlier than suggested 
by several pieces of evidence such as the observationally in- 
ferred SSFR-M S and F s -M s relation s at redshifts z < 1. 
As pointed out in lFirmani et al.l d2010l) and FA 10, the M s as- 
sembly of low-mass modeled galaxies follow, in first instance, 
the mass assembly of their corresponding halos, and in the 
ACDM scenario, low-mass halos assemble on average ear- 
lier than more massive halos (upsizing). The different astro- 
physical processes followed in current models and numerical 
simulations of galaxy formation and evolution (gas cooling 
and infall, SF, SF feedback, SN-driven galaxy outflows, etc.) 
makes the growth of M s deviates from that of the correspond- 
ing dark halo but not as strongly as observations apparently 
suggest (downsizing; see Fig. 4 in FA10 for a comparison 
of the average trends of halo MAHs and empirically inferred 
average galaxy M s tracks). 

If the empirical picture of sub-M* central galaxies assem- 
bly is confirmed (see §§4.2 for a discussion in current uncer- 
tainties), then the M s growth of central galaxies less massive 
than M s (z = 0)~ 3 x 10 10 M seems to require a signifi- 
cant delay with respect to the evolution of their correspond- 
ing ACDM halos; besi des, the smaller the g alaxy the longer 
such a delay (see a lso iNoeske et al!~ll2007bl) . According to 
iBouche et alj d2010l) . this mass-dependent delay in the M s as- 
sembly can be explained by a halo mass floor M m [ n ps 10 11 
Mq, below which the halo-driven gas accretion is quenched. 
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APPENDIX 

THE COMPILATION OF OBSERVATIONAL WORKS 

In Fig. [5] the simulation results in the SSFR-i\/ s plane are compared with a large body of observational inferences at four z 
bins compiled from the literature. Following, a brief description of these inferences at each z bin is presented . 

a) z = panel: The linear fit to the SSFR-M, relation and its intrinsic width as determined in lSalim et al.ld2007l) for the sub- 
sample of normal star-forming galaxies (optical emission lines are used for classifying galaxies) in a volume-corrected sample of 
around 50,000 SDSS galaxies are shown. The dust-corrected SFRs were obtained by fitting (GALEX) UV and SDSS photometry 
to a library of dust-attenuated SPS models. While the fraction of galaxies with negligible Ha emission-line detection in their 
sample is significant, most of them are luminous (massive) ones. Therefore, their SSFR-A/g relation for 'star-forming galaxies' 
at low masses (M s ^2x 10 10 M Q ) is almost the same if galaxies with no Ha detection are included (compare the mentioned 
fit with the one carried out by the same authors for the entire sample, where a Schechter function was used for the fitting, long- 
dash ed line; the characteristic mass is M = 10 1103 M Q ). Results similar to lSalim et al l (|2007|) were found by other authors 
(e.g JSchiminovich et ai1l2007h . 

b) z = 0.33 and z = 1.0 panels: The median a nd standard deviation of logSSFR vs logAf s of the sequence of star-forming 
field galaxies as reported in INoeske et al. I d2007bl) are plotted in these panels (crosses with vertical error bars). The data used 
by these authors consisted of 2905 galaxies out to z = 1.1 from the All- Wavelength Extended Groth Strip International Survey 
(AEGIS). The SFRs for galaxies with robust 24 /im detections were derived from 24 /im luminosity + non-extinction corrected 
emission lines, and for galaxies below the 24 /j,m detection, from extinction-corrected emission lines. The symbols show only 
the M s range where the sample is > 95% complete in the redshift bins 0.20 < z < 0.45 and 0.85 < z < 1.10, respectively. 
For masses smaller than the sample completeness (M s ~ 10 s M and ~ 10 9 M Q , respectively) the derived SSFRs continue to 
increase on average for lower values of M s . Around 30% of galaxies in their sample at all z have not robust 24 pan or emission- 
line detections, which implies very low SFRs. The inclusion of these galaxies in the SSFR-Tl/, relation, decreases the overall 
median SSFRs and it would be lower than the one shown in Fig. However, almost all of these galaxies (> 90%) are in the red 
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TABLE A 1 

Summary of sources and corrections 



Source 


SFR tracer 


Tft-TP 
llvlr 


ivi s onsei 


*vT7P nf'lVot' 

orK OIIScl 


Salim et al. (2007) 


Luv, H a 


Chabrier (2003) 


0.0 


0.0 


Bell et al. (2007) 


LjR + Ljjy 


Chabrier (2003) 


0.0 


0.0 


Noeske et al. (2007) 


L24fj,m + emission lines 


Kroupa (2001) 


0.0 


0.0 


Santini et al. (2009) 


Lir + L uv 


Salpeter(1955) 


-0.25 


-0.176 


Damen et al. (2009) 


Lir + L uv 


Kroupa (2001) 


0.0 


0.0 


Dunne et al. (2009) 


Ll400MHz 


Salpeter (1955) 


-0.25 


-0.176 


Kajisawa et al. (2010) 


Luv (corrected) 


Salpeter(1955) 


-0.25 


-0.176 


Karim et al. (2010) 




Chabrier (2003) 


0.0 


0.0 


Rodighiero et al. (2010) 


Lir + Luv 


Salpeter (1955) 


-0.25 


-0.176 


Gilbanketal. (2011) 


[OII]-(H a ) 


Baldry & Glazebrook (2003) 


-0.08 


-0.18 



"Shift in dex to convert M s to the Chabrier IMF 
* Shift in dex to convert SFR to the Chabrier IMF 

sequence and have early-type morphologies, being besides relatively massive (iNoeske et alJl2007atlNoeske et al. Il2007bl) . 

In these panels are als o plotted the avera ge data (from stacking) taken from the SSFR vs M s plots at 0.2 < z < 0.4 and 
0.8 < z < 1.0 given in iBell et al. I ( 120071) for 'all galaxies' (open stars) and for blue galaxies only (blue asterisks). Results 
correspond to the Chandra Deep Field South (CDFS) sample; the COMBO-17 survey in conjunction with Spitzer 24fim data 
were used for estimating M s and SFR (for the latter, UV and 24pm luminosities were used). In the low-mass side (M s < 3 10 10 
M s ) and for all z bins, blue galaxies dominate in such a way that the average SSFRs fo r the total sample ( blue + red galaxies) is 
close to the average of the blue cloud sample. The mass completeness limits reported in lBell et al. I(l2007l) are indicated with the 
dotted vertical arrows. 

c) z = 0.33, z = 1.00, and z = 1.50 panels: ISantini et al.l (120091) presented least-square linear fits to the SFR-M, relations 
(SFR=AM S ' 3 ) obtained from the GOOD-MUSIC catalog, which has multi-wavelength coverage from 0.3 to 24 pm (complete 
SED fittings and UV-24 /im luminosities were used for estimating SFRs). We plot in Fig. |5]these fits but divided by M s (blue 
lines) for the redshift intervals 0.3 < z < 0.6 (second panel), 0.6 < z < 1.0 and 1.0 < z < 1.5 (third panel), and 1.0 < z < 1.5 
and 1.5 < z < 2.5 (fourth panel). 

Similarly, the fits to the SSFR-i\/ s relations obtained in lRodighiero et ail (1201 Ol) for the redshift bins 0.3 < z < 0.6 (second 
panel), 0.6 < z < 1.0 and 1.0 < z < 1.5 (third panel), and 1.0 < z < 1.5 and 1.5 < z < 2.5 (fourth panel) are plotted (black 
solid lines). These authors used deep observations of the GOODS-N field taken with PACS on board of the Herschel satellite, 
which allows them to robustly derive the tot al infrared luminosi ty of galaxies, and combine w ith the ancillary UV luminosities 
in order to calculate the SFRs. In both cases ISantini et al.l (120091) and lRodighiero et aH (1201 Ol) . we plot their fits only for masses 
where samples are complete. 

The linear fits to radio-sfacfc/ng-based measurements of the SSFR as a function of M s recently presented bv lKarim et al.l(l201 ll) 
are also plotted in these three panels (long-dashed red line). The fits correspond to star-forming galaxies from a deep 3.6-/^m- 
selected sample of > 10 5 galaxies (0.2 < z < 3) in the 2 deg 2 COSMOS field, and are shown only for galaxies above the 
completeness limit at the given z. The fits to the entire sample are steeper than those to only star-forming galaxies, in such a way 
that at low masses both intersect, but at larger masses, the former lie below the latter. The image-stacking technique applied in 
order to increase the signal-to-noise ratio in the VLA 1 .4 GHz radio continuum observations allows one to estimate only average 
SFRs of the stacked population; they cannot shed light on the intrinsic dispersion of individual sources. 

Also shown in these panels are the average SSFRs corresponding to the lowest mass bins -where given samples are still 
complete- from: 

Damen et al. (2009b).- The mass bins of 0.5 dex width in each panel correspond to z sa 0.4, 1.0 and 1.4, respectively; 
the horizontal error bar indicates the width of the M s bin, 0.5 dex, and the vertical error bar represents the boots trapped 68% 
confid ence level on the average SSFR in the bin; the FIREWORKS catalog for the GOODS-CDFS generated by IWuvts et al.l 
was used to infer M s and SFRs (infrared and UV luminosities were combined to derive the latter and no selection on 
color/morphology was applied; i.e., all galaxies were included). 

Dunne et al. (2009).- The lowest complete mass bins of 1 dex width correspond to z = 0.50, 0.95, and 1.40, respectively, 
with the average SSFRs and their scatters shown with the error bars centered in the median value of the mass bin; stacked deep 
radio mosaics of if— band selected galaxies from the UKIDS S were us e d to d etermine M s and SFRs -rest-frame 1400-MHz 
luminosities were used to derive the latter; the case using the lBelletal](l2003T) conversion and applied to the whole sample is 
reproduced here. 

Kajisawa et al. (2010). -The lowest complete mass bins of 0.5 dex width correspond to the median dust-corrected SFRs, where 
the rest-frame 2800 A luminosity was used as indicator. This UV luminosity was estimated by the best-SED fitting of the multi- 
band photometry (UBVizJHK, 3.6 pm, 4.5 pxn, and 5.8 /im) with the BC03 models. We reproduce the reported values at redshift 
bins of 0.5 < z < 1.0, and 1.0 < z < 1.5. 

After the completion of this paper, appeared the work by iGilbank et al.l (1201 ll) . where the SSFR-M, relation for star-forming 
galaxies at z ~ 1 is presented down to masses smaller than those previously determined (their linear fit is plotted in panel 
z = 1 of Fig. [5] with a long-dashed purple line). For the low-mass side, they analyzed a sample of 199 low-mass galaxies 
(10 8 - 5 £ M S /M Q £ 10 9 3 ) from the field CDFS, with spectroscopy redshifts determined in the range 0.88 < z < 1.15 and a 
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sampling completeness > 80% dGilbank et al.ll2010h . Stellar masses were determined by SED-fitting of the photometry, at the 
given z, using a grid of PEGASE.2 models. The [Oil] line is used as the SFR indicator (actually the Kennicutt 1998 SFR-Ha 
relation is used by assuming an [OII]/Ha ratio of 0.5), empirically corrected as a function of M s in order to correct for extinction 
and other systematic effects. Their correction implies on average a lower SF R than previously det ermined, by using for example 
[Oil] (non-corrected) + 24/im (see for a discussion §§4.2. 1). The relation by iGilbank et al.l (1201 ll) is indeed the one that implies 
the lowest estimates for SSFR among all those compiled in Fig. [5]at zwl, These authors estimated also the local SSFR-M, of 
blue galaxies by usi ng SPSS Stripe 82 and the Ha line as the SFR indicator. The linear fit to their results lies below the average 
relation reported bv lSalim et all (I2007L by 0.25 dex at 10 9 M ). 

In Table I All a summary of all the sources used in Fig. |5]and described above is presented. In all the cases, h = 0.70 or 0.71 
was assumed. Columns 2 and 3 report the main SFR indicator(s) and IMF used in each source, r espectively. Col umns 4 and 
5 give the average shift s in dex applied h ere respectively to M s and SFR in order to correct for a IChabrierl (12003b IMF when 
necessary. According to lBell et aLl(l2007l) . the stellar masses and SFRs calculated with lChabrieri (120031) and lKroupal (1200 ll) IMFs 
(when using L jr + Luv as indicator) are consistent to within ^ 10%, so we do no apply corrections for those cases when a 
iKroupal J200l IMF was used (see for similar conclusions Kari m et al. 2010 a nd Salim et al. 20 07 for the L\^qqmhz an d H a 
indicators , respective l y). Fo r the corrections in cases when the (Salpeter] |1955l) IMF was used, see ISantini et al. ( 20 091). For the 
iBaldrv & Glazebroo k (2003) IMF, we have used the corrections given in IGilbank et alj d201 ll) for passing to a IKroupal ( 120011) 
IMF. We are aware that all these corrections are yet uncertain, and could actually vary from indicator to indicator and with z. In 
any case, they are relatively small, in particular for the SSFR. 
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